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Abstract 

We show that within MSSM singlet extensions the experimental hints beyond 
the standard model from the Fermi LAT telescope as well as from the LHC can be 
explained simultaneously while being consistent with all experimental constraints. 
In particular we present an example point which features a ~ 130 GeV lightest 
neutralino with an annihilation cross section into photons consistent with the 
indication from the Fermi satellite with simultaneously the right relic abundance, 
a continuum photon spectrum consistent with observation, direct detection cross 
section below the experimental limits, electroweak observables consistent with 
experiment and a 125 GeV light Higgs boson with a slightly enhanced h — > 77 
rate. 



1 Introduction 



2012 has been a very intriguing year regarding hints for new physics, both at the LHC 
and the Fermi large area telescope. While the data from ATLAS [I] and CMS |2) feature 
a new bosonic state with mass m ~ 125 GeV which is consistent with the expectation 
from a standard model (SM) Higgs, both experiments see indications for an excess in 
the diphoton channel. This potential enhancement in the diphoton rate has attracted 
much attention recently, see e.g. [3 20 . In contrast the diboson decays into WW* and 
ZZ* seem to be in accord with the SM expectation, which make explanations of the 
enhanced diphoton rate due to an increased partial decay width particularly appealing. 

A similarly exciting topic this year have been the hints for a 7-ray line in the Fermi 
LAT data as reported in 21,22 . 7 ray lines are considered the smoking gun of an- 



nihilating dark matter, as astrophysical processes able to induce line-like features are 
(see however 



23 241). Intriguingly, the morphology of the excess is consistent 



very rare 

with the expected distribution of dark matter up to a small offset from the galactic 
center j25,26|J] Currently the data are being re- analysed by the Fermi collaboration, 
and there seems to be an indication of a line-like feature at a slightly higher energy of 
135 GeV [28], where the shift results from a reprocessing of the data. The statistical 
significance of the excess found by the Fermi collaboration is, however, not as high as 
claimed in 



22,25,29 



although this also depends on the target region considered. A 
line feature has also appeared in the 7 ray data of the earth limb, raising some concerns 
about an instrumental effect 28,30,31 . Radio telescopes might help to confirm or rule 
out the dark matter interpretation of the line soon 32 . While the origin of the 7-ray 



line from the galactic center still has to be clarified, a noticeable amount of theoretical 
interest has been triggered 33 -49 . If interpreted in terms of dark matter, the 7-ray line 



requires a rather large annihilation cross section into photons, (av ) 77 ~ 10 



-27 3 - 

cm s 



if one assumes an Einasto profile 22 . Such a large cross section was found to be very 
difficult to accommodate in particle physics models, especially as the Fermi data are 
consistent with pure background at lower energies, i.e. competing annihilation channels 
must be sufficiently suppressed |42|[47|[50|[5l] . 

While one may tackle each of these experimental anomalies individually, it is in- 
triguing to speculate about a possible common origin. In this article we show that the 
signals from both Fermi and the LHC can be explained simultaneously within singlet 
extensions of the MSSM, while being consistent with all experimental constraints. Sin- 
glet extensions are particularly interesting given the observed value for the Higgs mass 



Recent numerical simulations indicate that such a small offset can indeed be realized in realistic 



models of galaxy formation 27 
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because the electroweak fine-tuning can be substantially alleviated in these models. A 
somewhat generalised version of the NMSSM, the GNMSSM [52], which is based on a 



discrete R symmetry 53 , 54 , was found to be particularly promising in this context 55 



In 19 it was found that in this setup the coupling of the CP even neutral light Higgs to 



light charginos can be strongly enhanced, leading to a sizeable increase in the h — > 77 
rate. Interestingly the same coupling can lead to a rather large neutralino annihilation 
cross section into photons as indicated by the Fermi data, while being compatible with 
bounds from direct detection, electroweak precision observables, the continuum photon 
spectrum and with a 125 GeV Higgs with an enhanced diphoton rate. 

This article is organised as follows: In the next section we will briefly review some 
aspects of the GNMSSM. In section[3]we will then discuss neutralino annihilation within 
this framework, before we come to constraints arising from the requirement of the correct 
relic abundance and the continuum photon spectrum in section |4} Constraints from 
direct and indirect detection experiments are analysed in section [5j Section [6] is then 
devoted to a thorough numerical study of a benchmark point, while section [7] contains 
our summary. Some useful information about the GNMSSM is collected in the appendix. 



2 The GNMSSM 

As a framework, we consider the GNMSSM, a generalised version of the NMSSM, which 
has a superpotential of the form 

W = W Yukawa + ^S 3 + {fj, + \S)H u H d + l -^ s S 2 . (1) 

Here WVukawa are the MSSM superpotential terms generating the usual Yukawa couplings 
and we used the freedom to shift the singlet S to set a potential linear term in S to zero. 
This superpotential has additional explicit mass terms n and fi s which are not present 
in the Z 3 symmetric NMSSM which is usually considered (for reviews of the NMSSM 



see e.g. 56,57|). While the apparent un-naturalness of these additional mass terms has 
prevented a larger community from studying the phenomenology of the GNMSSM, it 
has recently been realised that exactly this structure naturally arises from an underlying 



R symmetry as discussed in 54 . The fact that this R symmetry also eliminates the 
dangerous dimension four and five baryon- and lepton-number violating terms and avoids 
destabilising tadpoles and domain wall problems makes it a more promising starting 
point than the symmetric NMSSM. 

The soft SUSY breaking terms associated with the extended Higgs sector of the 
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GNMSSM are given by 



soft 



Z\ 12 , 2 1112, Z II I 

™>l + m hjhu\ +m hd \h d \ 



bfi h u h d + XA x sh u h d + ^kA k s 3 + ^b s s 2 + £ s s + h.c. 



(2) 



The resulting mass matrices as well as the relevant couplings for our discussion are given 
in appendix [A] and [B} The field content is the same as in the NMSSM. In comparison 
with the MSSM there is an additional singlet fermion which mixes into the neutralino 
sector and an additional complex scalar which mixes into the Higgs sector. For more 



details on the model, see 52 , 55 
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Figure 1: In the GNMSSM the A-coupling drives the annihilation of neutralinos into pho- 
tons (left), enhances the Higgs decay rate into photons (middle) and provides an additional 
contribution to the Higgs boson mass (right). 



The Higgs sector of the GNMSSM has been discussed in 19152 55 . As in the NMSSM 



there is an additional tree-level contribution to the lightest Higgs mass, which is large for 
small tan (3 and large A, allowing to evade the upper bound of rrih < Mz present in the 
MSSM. The fact that radiative corrections due to top/stop loops are no longer needed 



drastically reduces the electroweak fine-tuning in MSSM singlet extensions 58 -62 for 
large values of A. Given the Higgs mass of 125 GeV the GNMSSM has been shown 
to be particularly interesting in this context 52 , 55 . One additional advantage of the 



GNMSSM is the fact that for small v s no tuning in the Higgs mass matrix is required 
to avoid large doublet-singlet mixing^} It has also been argued that allowing even larger 
values of A, relaxing the condition that it remain perturbative up to the GUT scale, 
leads to an additional reduction in the fine-tuning 



63 



2 In the NMSSM it is common practice to tune the parameter A\ such that the off-diagonal entry 
in the mass matrix is close to zero, avoiding doublet-singlet mixing and allowing for a correspondingly 
larger Higgs mass. 
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As discussed in [19] scenarios with large A may also accommodate an enhanced 
diphoton rate as observed by ATLAS and CMS: in addition to the dominant W and top 
loops present in the SM, new sizeable contributions from charged Higgs and chargino 
loops can arise, if these states are light. While a light charged Higgs is generically 
challenged by the measurement of b — > 57, a light chargino is perfectly viable and an 
interesting possibility. In the next section we will show that large A and light charginos 
can also lead to a large annihilation cross section into photons, see also figure [TJ 



3 Neutralino annihilation into photons 



The mass matrix as well as the relevant couplings of the GNMSSM neutralino sector can 
be found in the appendix. Compared to the MSSM, the GNMSSM offers two alterna- 





Figure 2: Relevant Feynman diagrams for annihilation of the lightest neutralino into two 
photons within the GNMSSM. 

tive possibilities to achieve a large annihilation cross section (av ) 77 . The corresponding 
Feynman diagrams are shown in figure [2J In case Xi carries a sizeable singlino frac- 
tion, the chargino/ charged Higgs loop on the left experiences a drastic enhancement 
compared to the MSSM through the A-coupling. However, even for A > 1, a cross 
section large enough to explain the Fermi line requires not only the charginos but also 
the charged Higgs bosons to be very light. This typically causes problems with flavour 
observables, in particular b — > 37. Therefore, we will concentrate on the second pos- 
sibility in the following, shown in the right panel of figure [2j The diagram contains a 
pseudoscalar Higgs Aj in the s-channel, i.e. the cross section can be enhanced in the 
vicinity of the pseudoscalar resonance (see also [44] ) . 



Analytically, one finds (see e.g. 64, 65] 
Mrr= ~T^t * 



(3) 
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Here G^v /!! = 9~o~a A — Q-o-o A and 9o±^±4 

i'XiXi^l ^XiXi^l ^XiXi^l X, ^1 

of the trilinear couplings g^o^o Al and g^~ ±Ai is given in appendix 



The general form 
B Throughout this 

article we will assume A and k to be real. 

An interesting limit arises if we assume that Xi and A.\ are dominantly singlet-like. 
In this case, we obtain g^a^o Al = y/2n. Further, only the higgsino-like chargino couples 
to A\ with g~±~± A = A/a/2 for the lighter chargino being a pure higgsino. Setting as 

X.\ A,\ "1 

an example m^± = 1.5 m%o, we obtain in this limit 



(cru) 77 ~ (6 ■ 10 



" 28 cm 3 s" 



l ) ■ AV 



(100 GeVf 
4 m 2 - m A 

Al 



Al 



130 GeW 



(5) 



This shows that a cross section large enough to explain the Fermi line can indeed be 
realised if we allow for a mild tuning of rriA 1 - In the above example, m Al has to be 
within the range rriA 1 — 240 — 280 GeV (if we assume A, ft < 1). The tuning is, how- 
ever, substantially less severe than in the NMSSM case Ell. The reason is that in 
the NMSSM it is not possible to obtain a dominantly singlet like pseudoscalar and a 
singlino like LSP with the desired masses simultaneously. Instead the lightest neutralino 
has to be predominantly bino like, leading to much smaller effective couplings. Large 
enough annihilation cross sections can only be achieved very close to the pseudoscalar 
resonance, which in turn requires substantial tuning, in particular, when the most recent 



XENON100 bounds 66 are taken into account, which appeared after 44 was published. 



In addition to the 7 ray line from annihilation into two photons, there is a second 
line from the process X1X1 ~~ %1 at a slightly lower energy 



m x1 I 1 



Ml 

x?. 



(6) 



As has been shown in 25 such a double line structure even slightly improves the fit to 



the Fermi data compared to a single line. The cross section (<7v)z 1 can be calculated with 
the formulas presented in 67 . We find that for a pure higgsino (av ).z 7 /(cn;) 77 ~ 0.6. 
The wino, however, has a stronger coupling to Z bosons than the higgsino. Therefore, 
chargino mixing tends to increase the relative importance of the ^-channel. 
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4 Continuum photons and relic density 



Even if dark matter annihilations induce a 7 ray line consistent with Fermi, it remains 
quite challenging not to overproduce continuum 7s by competing annihilation processes. 
Irrespective of whether dark matter is produced thermally or non-thermally, its present- 



day annihilation fraction into 7s must satisfy 47 



Br = > 1(T 2 (7) 

tfr 77 ^ „iu . {() 

For the winos and higgsinos of the MSSM, this fraction is in the range Br 77 = O(10~ 3 ). 
Therefore, the low energy 7 data disfavour an explanation of the Fermi line within the 



framework of the MSSM 47 ,50 



In the GNMSSM there is a simple possibility to suppress the continuum 7s. As 
in 44 , we take the lighter pseudoscalar A\ to be an almost pure singlet. With the 
appropriate choice of 6/i and b s it is simple to arrive at a situation where the singlet 
pseudoscalar remains light while the MSSM pseudoscalar becomes heavy and decouples. 
In this case, any tree level annihilation process into quarks through an intermediate 
Ai is suppressed by the mixing angle between singlet and MSSM pseudoscalar. More 
specifically, for A ~ 1 the effective coupling of Ai to photons becomes comparable to the 
tree-level coupling to bottom quarks for a doublet fraction at the level of 0.1%. This 
can easily be achieved for a sufficiently heavy A 2 without requiring any cancellations in 
the pseudoscalar mass matrix. 

Competing annihilation processes into SM states which do not proceed through A\ 
can be suppressed if Xi is dominantly singlino-like. A singlino-like LSP arises if the 
gauginos and higgsinos are sufficiently heavier than the singlino. With a singlino-like x? 
and a singlet-like Ai one naturally obtains Br 77 ^> 10~ 2 . 

Within the GNMSSM it is also possible to realise thermal dark matter and the 7 
ray line simultaneously. Thermal production requires a total neutralino annihilation 
cross section (av )fo ~ 2 ■ 10 -26 cm 3 s _1 at the time of freeze-out. This can e.g. be 
achieved through a subdominant wino admixture in xi which induces annihilation into 
W bosons. Alternatively, a thermal cross section can be realised by the annihilation of 
Xi into bottom quarks or gluons through the MSSM admixture of A\. 

In figure [3] we present the 7 ray flux for a specific benchmark choice of the GNMSSM 
parameters which can be found in table [TJ The benchmark scenario will be discussed in 
more detail in section [6] For the dark matter density distribution we assume an Einasto 
profile with the parameters as given in [22]. The Fermi data are taken from region 4 
(source) of the same reference. The 7 ray lines from dark matter annihilation into 77 
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Figure 3: 7 ray flux for the benchmark scenario of table [TJ The Fermi data are taken from 
(region 4, source). The dark matter induced flux comprises the lines from annihilation into 77 
and jZ as well as the continuum photons from the fragmentation and decay of the accompa- 
nying annihilation products. The background is modelled with a power law. 



and Z7 have been convoluted with the Fermi energy resolution as extracted from table 3 
in 



22 . The 7 ray background is modeled with a featureless power law with index —2.5. 
It can be seen that the double line at E — 130 GeV and E = 114 GeV (from annihi- 
lation into 77 and Zj respectively) gives a very good fit to the data. The continuum 7s 
which - in the benchmark scenario - mainly arise from the jyjy-channel are sufficiently 
suppressed, i.e. they are nicely consistent with the low energy data. 



5 Direct and indirect detection constraints 

Let us now briefly discuss constraints on the model which arise from direct dark matter 
searches. The spin- independent cross section of the lightest neutralino with nucleons 
o"^ 1 is typically dominated by exchange of the light scalar Higgs h\. It can be written in 
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where f2, and Jq denote the up-, down-, strange-quark and gluon contributions 
to the nucleon mass m n which we take from 68 . For simplicity, we have applied the 
decoupling limit on the MSSM Higgs fields such that the effective neutralino quark 
coupling divided by the quark mass f q is universal among the quark families. The latter 
takes the form (see e.g. [69]) 



\f q 



Gi 



m 



hi 



1 



7 H* 

'13 



(9) 



where Gf is the Fermi constant and the singlet fraction of the light Higgs hi (see 
appendix A). The neutralino Higgs coupling is given as g^o^ = g? ^ ()h ^ + gjp~o hi with 
the left ana right couplings being defined in appendix [B] Note that they depend strongly 
on the composition of the lightest neutralino. 

The relevant limit on the cross section a^ 1 is set by the XENON100 experiment, 

130 GeV corresponds to 



which for a neutralino mass m^o 

A l 



66 



a 



si 



< 3 • 10" 45 



cm 



(10) 



By use of Q we can translate this into a limit on the coupling g-g~ , . Assuming that 

hi is dominantly SM like, we estimate |<7£?x?/iil ^ O.OSj 4 

As described in section |4j we are mainly interested in the case where Xi is singlino-like 
such that the production of continuum 7s is suppressed. However, in order to enhance 
annihilation of Xi m to photon pairs, it is favourable to have light (charged) higgsinos. 
Thus Xi always contains a non-negligible higgsino admixture. In this case ^°x?fti typi- 
cally receives comparable contributions through the A as well as the k coupling bearing 
also the possibility of (partial) cancellations. Still, we find that for A ~ 1 the higgsino 
fraction of Xi should not exceed ~ 10% in order to satisfy the XENON bound. 

Further constraints on the model arise from the neutrino searches by Super-Kamio- 
kande 



70,71 



and IceCube [72] which aim to detect the annihilation of dark matter 
particles in the sun. They provide especially strong bounds on the spin-dependent cross 
section of WIMPs with protons. Assuming m^o ~ 130 GeV, the relevant upper limit 

if the LSPs dominantly annihilate into W bosons. The 



reads <Tp D 



4 • 10 



-40 cm 2 



70 



3 We neglect the small differences between proton and neutron. 

4 Note, however, that the uncertainties in the nucleon composition may affect this constraint. 
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constraint, however, gets significantly weaker for annihilation channels which induce a 
softer neutrino spectrum. In the considered region of parameter space, the leading spin- 
dependent WIMP proton cross section <7p D arises from Z exchange. It scales with the 
higgsino components of Xi, more specifically <rjp oc | ^V 13 1 2 — | A^"i 4 1 2 (cf. appendix p 



While this cross-section can be sizeable, there generically occur cancellations due to 



\N- 



13 1 



\N- 



14 • 



The most dangerous situation arises if Xi is strongly mixed among the 
different states as in this case iV" 13 and iV" 14 typically get split. Then, a crp close to the 
current experimental bounds may be generated. 

We find that the remaining constraints from indirect dark matter detection, arising 
e.g. from 7 ray searches in the Milky Way satellite galaxies by Fermi 73 or antiproton 



searches by PAMELA 74,75 and BESS-Polar II 76,77 , are in general weaker than 



those from the continuum 7s studied previously. To illustrate that in the GNMSSM, 
the Fermi line can be explained while all direct and indirect detection constraints are 
satisfied, we provide an explicit example in the next section. 



6 A benchmark scenario 



After we have gained some analytical understanding, we now turn to a full-fledged nu- 
merical analysis. For this purpose we use the SPheno version 78,79 for the GNMSSM 
created by SARAH 



80 83 which has been presented in 55 



complete one-loop calculation of all SUSY and Higgs masses 84 



dominant two-loop corrections for the scalar Higgs masses 86 - 89 



This version performs a 
and includes the 
In addition, it cal- 



85 



culates the decay widths and branching ratios of all SUSY and Higgs particles. In the 
Higgs sector the decays are calculated with the following precision: the channels with 
two SUSY particles, SM leptons or SM vector bosons in the final state are calculated at 
tree level. In contrast, for quark final states the dominant one-loop QCD corrections due 
to gluons are included 



90 



For the decays into two photons and two gluons induced at 
one-loop level all possible leading order contributions are included. In addition, for the 
CP even Higgs also the dominant NLO QCD corrections are added 90 . Furthermore, 



this SPheno version also includes routines to calculate b — > 57, 6p and g — 2 which have 
been used to check possible constraints from these observables. These calculations are 
performed in the GNMSSM with the same precision as described in Ref. 79 for the 
MSSM. 

For the calculation of the relic density of the lightest neutralino as well as to obtain 
(cru) 77 we have used MicrOmegas 91 93 . For this purpose we created model files for 
CalcHep 94 with SARAH. These model files include optionally also the effective interac- 
tions /ij77 and ^77. The numerical values for these operators as well as of all other 
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Figure 4: Left: The present day diphoton rate (cu} 77 (solid line) and the total cross section 
{(tv) (dashed line) as function of the light pseudoscalar mass m^. Right: dependence of the 
dark matter relic density VLh? on the pseudoscalar mass. The other parameters are those of 
the point given in table [T} The dotted, vertical line indicates the pseudoscalar mass for our 
benchmark point. 



parameters are read from the spectrum file written by SPheno using the SLHA+ function- 



ality of CalcHep 95 . We also used MicrOmegas to calculate the continuous 7 spectrum 
which has already been discussed in section |4j 

In table [T] we show a benchmark point with all the desired features: the light Higgs 
mass is close to 125 GeV and the branching ratios into two photons is enhanced by 20% 
because of the chargino loop contributions. The LSP is the lightest neutralino which is 
mostly singlino-like and has a mass of 130 GeV. A mostly singlino LSP can be achieved 
by appropriate values for M 1; M 2 and /i e g. While M\ plays only a subleading role as long 
as the bino is heavier than the singlino, the choice of // e g has to be done more carefully: 
on the one hand it should not be too small in order to suppress the mixing between the 
singlino and higgsino, because a sizeable higgsino fraction is often in conflict with direct 
detection measurements. On the other hand, light charginos with a sizeable higgsino 
fraction are needed in order to enhance the loop contributions to h — > 77 and (en;) 77 . 
For the benchmark point a light chargino has been realized by a comparably small value 
of Mii which however leads to a large mixing in the chargino sector. Note that the 
benchmark scenario is consistent with direct chargino searches at the LHC as the latter 
have only gained sensitivity to spectra with m^o < 100 GeV 97 [98] Pj The correct relic 
density is obtained mainly via annihilation into W + W~ , while today s annihilation into 
photons mainly proceeds via the pseudoscalar exchange. This pseudoscalar is nearly a 

5 The sensitivity of chargino searches increases if there exists a slepton with mass between m~+ and 
m-^o . This is, however, not the case in our benchmark scenario. 
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pure singlet and it is not necessary to be very close to the resonance: even with a mass 
more than 10 GeV away from the resonance the diphoton rate is enhanced to a level 
sufficient to explain the tentative Fermi line. This is a big improvement in comparison to 
the NMSSM where one has to be usually very close to the resonance: for cases with the 
pseudoscalar component of the singlet in the correct mass range, the singlino fraction of 
the LSP is very small and the coupling between both is highly suppressed. This is not 
only a drawback of the NMSSM with respect to the needed fine-tuning, but, what is even 
more important, this scenario is also under big pressure from direct detections bounds. 
The dependence of (av) on is shown in figure [4] (left), with all other parameters 
as given in table [TJ In the right panel we also show the relic density Vlh 2 , which in the 
region of interest is hardly effected by the pseudoscalar resonance. 

In order to illustrate the constraints on the parameter space of the GNMSSM which 
arise from direct and indirect detection, we show the results of a scan over /i and M 2 in 
figure |5} For each combination (/z, M 2 ) we adjusted the parameters fi s and b s such that 
m^o and rriA 1 remain fixed at their values from table [TJ The remaining parameters were 
taken from the same table. 



260 FT" 




250 260 270 280 290 300 310 320 
fi [GeV] 



Figure 5: Parameter Scan in the GNMSSM (see text) with contours referring to the value of 
(<tv) 77 . The benchmark scenario of table [TJ is indicated by the blue dot. In the gray band, 
the thermal relic density of the lightest neutralino is consistent with the observed dark matter 
density. The yellow region is exluded by the XENON100 direct dark matter search, the red 
region is excluded by the Super-Kamiokande limit on the spin-dependent cross section of Xi- 

It can be seen that, while Xi is dominated by its singlino component in the depicted 
parameter space, a rather low M 2 is favorable in order to satisfy the direct detection 
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constraints. This is because the Higgs, which mediates the spin-independent interaction 
of Xi with nucleons, mainly couples to the higgsino and singlino components of Xi due 
to the relatively large A and k. An additional wino admixture tends to reduce cx^ 1 . At 
the same time, the wino component of Xi drives its annihilation into W pairs allowing 
for a thermal neutralino abundance which agrees with the dark matter abundance (gray 
band in figure [5J. On the other hand, a light wino induces a splitting between the 
higgsino components of the lightest neutralino, N 13 and AT 14 . This, in turn, enhances 
the spin-dependent cross section <r^ D (see section 5]) resulting in the Super-Kamiokande 
excluded region at low and /ijj The annihilation of xi m t° 7-pairs is mediated by 
higgsinos in the loop, while winos do virtually not contribute due to the singlet nature 
of Ax. Therefore, (cn;) 77 grows with decreasing higgsino mass. Note, finally, that the 
direct and indirect detection constraints get weaker for smaller couplings A and k. In 
this case, however, one would need closer to the resonance in order to keep (av ) 77 
large. 



7 Summary 

In this article we have shown that within the GNMSSM, a generalised version of the 
NMSSM, the experimental hints beyond the SM from the Fermi LAT telescope as well 
as from the LHC can be explained simultaneously while being consistent with all experi- 
mental constraints. As in the NMSSM the superpotential term A SH u Hd plays a crucial 
role in this scheme: it drives the annihilation of the lightest neutralino into photons, 
induces new tree-level contributions to the mass of the light Higgs boson h and enhances 
the partial width of the decay h — > 77. To obtain a large enough annihilation cross sec- 
tion into photons, a very mild tuning of the pseudoscalar Higgs mass is required in the 
GNMSSM: it is sufficient to be within ~ 20 GeV from the resonance. This is in contrast 
to the NMSSM, where the tuning is very substantial, because unlike in the GNMSSM 
a mainly singlino-like LSP is not possible, leading to much smaller effective couplings. 
Hence while the coupling A is also present in the NMSSM, it is the additional flexibility 
in the mass spectrum of the GNMSSM which allows to simultaneously explain the Fermi 
7 ray line, ameliorate the little hierarchy problem, and explain a moderate excess in the 
77-channel which is indicated by the Higgs searches of CMS and ATLAS. 

To substantiate our claim we performed a thorough numerical analysis and presented 
an example point which features a ~ 130 GeV lightest neutralino with an annihilation 
cross section into photons consistent with the indication from the Fermi satellite with 

6 To obtain the Super-Kamiokande limit on the spin-dependent cross section, we have weighted o^ 1 
by the fractional annihilation of \i hito channels which induce hard neutrinos. 
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simultaneously the right relic abundance, a continuum photon spectrum consistent with 
observation, direct detection cross section below the experimental limits, electroweak 
observables consistent with experiment and a 125 GeV light Higgs boson with a slightly 
enhanced h — > 77 rate. 
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A Mass matrices 

If we decompose the complex Higgs fields and singlet after EWSB as 

H° = -j=(<t> d + v d + ia d ) (11) 
H° u = ^=( ( j )u + v u + ta u ) (12) 
S = -j= (0 S + v s + ia s ) (13) 

the mass matrices in the neutral Higgs sector read: 
• Scalar Higgs. Basis: {<fi d , M , <p s ) 

m dd = ^(4(2/i + y/2v 8 \)n* + 4(^/1 + (v 2 s + ^)a)a* 

+ 8m^+(^ + ^ 2 2 )(3^-^)) (14) 

m ud = -^-2v 2 s ^{\K^-A^(b^ +4^ U |A| 2 

- (gl + gi)v d v u - V2v s (2$t(\^ +2K(AA A ))) (15) 

m uu = - (a(2(jl + y/2v a \jfJL* + A(y/2v a fJL + (v 2 d + w 2 )a)a* + 8m 2 Hu 

-(£ + gl)(-K + ti)) (16) 

v u (Av s U^*) + V2(2M(XA X ^ +/i 5 A*) + v^A/i^) 



1 

m sd = --. 
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+ (-^ u // + ^ u a)a* 
™ ss = ^ (6^V| 2 + ^|A| 2 + ^|A| 2 + 2 (m| + |/i 5 | 2 + ^(6, 

Pseudoscalar Higgs. Basis: (a d: <J U , a s ) 

m dd = ^ (a(2/i + yfivsXjp* +4(V2v s ii + (y 2 s + vty a) A* 

+ 8<+(5i+3 2 2 )(-^ + ^)) 
m ud = ^m(bf?j +v s (2v s ^\k*^ + ^(29^/4) +23?(aA; 

™« u = ^ (4(2/i + v^a)// +4(^2^// + (f d + ^) A ) A * 
+ 8m k- (^i+^)(-^ + ^ 



m sd = ~v u (±v&(\iC) + V2( - 2K(AA A ) + /i 5 A*) + v^A/^) 
m su = -^(^(Ak*) + - 2ft (AA A ) + /i 5 A*) + v^A/^) 
m ss = \ (2v\\^ + v\\\f + ^|A| 2 + 2( - »(&*) + m| 

+ |^ 5 | 2 ) + 2v d vJt(\K^ + v / 2^(2$R(^) - 2&(«4 S ))) 

After eliminating the Goldstone mode, the mass matrix can be written as 
M 1;1 = (v s (V2(A x + pL S ) + v s k)\ + 2bfi) /sin(2/3) 

M 1>2 = ^(A A - \I~2~V S K - fJL a ) 

3 

M 2 ,2 = -26 s ^A^/t - (v 2 Xfi)/(y/2v s ) - v s K,fi s /V2 - V2£ s /v s 

V2 

14 



+ ((A\ + Us) Xv 2 cos (3 sin (3) /(V2v s ) + 2v 2 kX cos /3 sin (3 



(30) 



The mass matrices of the charged Higgs as well as of the neutralinos and charginos are 
given by 

• Charged Higgs. Basis: (H d ,H+'*) 



m dd = - 4 



v s X j /I* + Av s (y/2fj, + v s X^j A* + 8m 2 Hd 
+ 9l(-v 2 u + v 2 )+ g 2 2 {v 2 d + vi)) 
m ud = ^(2v s (y2(\^ s + \A x ^j + v s Xk*^ + 4bfi - v d v u (2\\\ 2 - gfj^j 
m uu = - (aV2v s \ii* + Av s (V2fi + v s X^j X* + 8m 2 Hu 
+ 8\fi\ 2 + (gl + gt)v 2 u - g\v\ + g 2 2 v 2 d ) 

Neutralino. Basis: W°, H° d , 



(31) 
(32) 

(33) 



/ M 1 




m N 





M 2 



-\g\Vd \g2v d 
\giv u -\g2Vu 



-\g\v d 
\giv d 




Chargino. Basis: (w~,H^ , (w + ,H+^ 



\g\v u 
-\g2Vu 
■^VsX-v 



"75^ 








(34) 



M 9 



m Ch 



V2 



g2Vv 



T2 g2v d T2 v s X + fi 



(35) 
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B Couplings 

g\, m =- ( - 92n: 2 n; 3 z^ - V^xn^z* - V2\n* 4 n; 5 z^ - g^N^Z* 

+ g 2 N; 4 N; 2 z£ 2 - V2\N* 5 N* 3 Z£ 2 + g2N* 2 N* 4 z£ 2 

+ 9i n; 1 (n; 3 z* - n; 4 z£) - V2\n; 4 n; 3 z£ 3 + 2V2 K n; 5 n; 5 z£ 3 

+ TV* [g^Zti - 92N; 2 Z^ - V2\ [n^Z* + N* 5 z£) ) ) 
9mo A H--(g^ m r\(i^j) (36) 



xtXj A h k 



g *tx-At =^ ( - 92U* n V* 2 Z£ 2 + E$ ( - g 2 V*Z^ + \V*Z£ 3 ) ) 



4^ = - ( %W* l( ^ (37) 



X°X°jh k 

9m° hk =\ ( - 92N* 2 n; 3 z» + V2\N* 5 N* 4 Z» + V2\N* 4 N* 5 Z» - g^N^Z* 

+ g 2 N; 4 N; 2 z» + V2\n* 5 n; 3 z? 2 + 92 n: 2 n; 4 z? 2 

+ g 1 N* 1 (N* 3 Z» - iV* 4 Zf 2 ) + V2\N* 4 N* 3 Zg - 2^nN^Zg 
9m° hk =-(9^ hk T\(i^j) (38) 



• x°x°z 



9 mz =(s^zY (39) 



References 

[1] ATLAS Collaboration, ATLAS-CONF-2012-093. 



16 



[2] CMS Collaboration, CMS-PAS-HIG-12-020. 

[3] U. Ellwanger, JHEP 1203 (2012), 044, [1112.3548]. 

[4] J.- J. Cao, Z.-X. Heng, J. M. Yang, Y.-M. Zhang, and J.-Y. Zhu, JHEP 1203 (2012), 
086, [1202.5821]. 

[5] V. Barger, M. Ishida, and W.-Y. Keung, (2012), 1207.0779. 

[6] M. Carena, I. Low, and C. E. Wagner, (2012), 1206.1082. 

[7] D. S. Alves, P. J. Fox, and N. J. Weiner, (2012), 1207.5499. 

[8] N. Bonne and G. Moreau, (2012), 1206.3360. 

[9] B. Bellazzini, C. Petersson, and R. Torre, (2012), 1207.0803. 
[10] M. R. Buckley and D. Hooper, (2012), 1207.1445. 
[11] H. An, T. Liu, and L.-T. Wang, (2012), 1207.2473. 
[12] A. G. Cohen and M. Schmaltz, (2012), 1207.3495. 

[13] A. Alves, A. Dias, E. R. Barreto, C. S. Pires, F. S. Queiroz, et al., (2012), 1207.3699. 

[14] A. Joglekar, P. Schwaller, and C. E. Wagner, (2012), 1207.4235. 

[15] N. Haba, K. Kaneta, Y. Mimura, and R. Takahashi, (2012), 1207.5102. 

[16] L. G. Almeida, E. Bertuzzo, P. A. Machado, and R. Z. Funchal, (2012), 1207.5254. 

[17] R. Benbrik, M. G. Bock, S. Heinemeyer, O. Stal, G. Weiglein, et al., (2012), 
1207.1096. 

[18] N. Arkani-Hamed, K. Blum, R. T. D'Agnolo, and J. Fan, (2012), 1207.4482. 
[19] K. Schmidt-Hoberg and F. Staub, (2012), 1208.1683. 
[20] M. Chala, (2012), 1210.6208. 

[21] T. Bringmann, X. Huang, A. Ibarra, S. Vogl, and C. Weniger, JCAP 1207 (2012), 
054, [1203.1312]. 

[22] C. Weniger, (2012), 1204.2797. 

17 



[23] S. Profumo and T. Linden, JCAP 1207 (2012), Oil, [1204.6047]. 
[24] F. Aharonian, D. Khangulyan, and D. Malyshev, (2012), 1207.0458. 
[25] M. Su and D. P. Finkbeiner, (2012), 1206.1616. 
[26] T. Bringmann and C. Weniger, (2012), 1208.5481. 

[27] M. Kuhlen, J. Guedes, A. Pillepich, P. Madau, and L. Mayer, (2012), 1208.4844. 

[28] A. Albert, Talk at the Fermi Symposium 2012. 

[29] E. Tempel, A. Hektor, and M. Raidal, (2012), 1205.1045. 

[30] A. Hektor, M. Raidal, and E. Tempel, (2012), 1209.4548. 

[31] D. P. Finkbeiner, M. Su, and C. Weniger, (2012), 1209.4562. 

[32] R. Laha, K. C. Y. Ng, B. Dasgupta, and S. Horiuchi, (2012), 1208.5488. 

[33] G. Chalons, (2012), 1204.4591. 

[34] G. Chalons and A. Semenov, JHEP 1112 (2011), 055, [1110.2064]. 

[35] A. Ibarra, S. Lopez Gehler, and M. Pato, (2012), 1205.0007. 

[36] E. Dudas, Y. Mambrini, S. Pokorski, and A. Romagnoni, (2012), 1205.1520. 

[37] J. M. Cline, (2012), 1205.2688. 

[38] K.-Y. Choi and O. Seto, (2012), 1205.3276. 

[39] B. Kyae and J.-C. Park, (2012), 1205.4151. 

[40] H. M. Lee, M. Park, and W.-I. Park, (2012), 1205.4675. 

[41] A. Rajaraman, T. M. Tait, and D. Whiteson, (2012), 1205.4723. 

[42] M. R. Buckley and D. Hooper, (2012), 1205.6811. 

[43] X. Chu, T. Hambye, T. Scarna, and M. H. Tytgat, (2012), 1206.2279. 

[44] D. Das, U. Ellwanger, and P. Mitropoulos, JCAP 1208 (2012), 003, [1206.2639]. 

[45] Z. Kang, T. Li, J. Li, and Y. Liu, (2012), 1206.2863. 

18 



[46] N. Weiner and I. Yavin, (2012), 1206.2910. 

[47] W. Buchmuller and M. Garny, (2012), 1206.7056. 

[48] M. T. Frandsen, U. Haisch, F. Kahlhoefer, P. Mertsch, and K. Schmidt-Hoberg, 
JCAP 1210 (2012), 033, [1207.3971]. 

[49] F. D'Eramo, M. McCullough, and J. Thaler, (2012), 1210.7817. 

[50] T. Cohen, M. Lisanti, T. R. Slatyer, and J. G. Wacker, JHEP 1210 (2012), 134, 
[1207.0800]. 

[51] I. Cholis, M. Tavakoli, and P. Ullio, Phys.Rev. D86 (2012), 083525, [1207.1468]. 

[52] G. G. Ross and K. Schmidt-Hoberg, (2011), 1108.1284. 

[53] H. M. Lee et al, Phys. Lett. B694 (2011), 491-495, [1009.0905]. 

[54] H. M. Lee et al, Nucl. Phys. B850 (2011), 1-30, [1102.3595]. 

[55] G. G. Ross, K. Schmidt-Hoberg, and F. Staub, (2012), 1205.1509. 

[56] M. Maniatis, Int. J.Mod.Phys. A25 (2010), 3505-3602, [0906.0777]. 

[57] U. Ellwanger, C. Hugonie, and A. M. Teixeira, Phys. Rept. 496 (2010), 1-77, 
[0910.1785]. 

[58] M. Bastero-Gil, C. Hugonie, S. King, D. Roy, and S. Vempati, Phys.Lett. B489 
(2000), 359-366, [hep-ph/0006198]. 

[59] R. Dermisek and J. F. Gunion, Phys.Rev. D73 (2006), 111701, [hep-ph/0510322]. 

[60] R. Dermisek, J. F. Gunion, and B. McElrath, Phys.Rev. D76 (2007), 051105, [hep- 
ph/0612031]. 

[61] R. Dermisek and J. F. Gunion, Phys.Rev. D76 (2007), 095006, [0705.4387]. 

[62] U. Ellwanger, G. Espitalier-Noel, and C. Hugonie, (2011), 1107.2472. 

[63] L. J. Hall, D. Pinner, and J. T. Ruderman, (2011), 1112.2703. 

[64] L. Bergstrom and P. Ullio, Nucl.Phys. B504 (1997), 27-44, [hep-ph/9706232]. 



19 



[65] Z. Bern, P. Gondolo, and M. Perelstein, Phys.Lett. B411 (1997), 86-96, [hep- 
ph/9706538]. 

[66] XENON100 Collaboration, E. Aprile et al., Phys.Rev.Lett. 109 (2012), 181301, 
[1207.5988]. 

[67] P. Ullio and L. Bergstrom, Phys.Rev. D57 (1998), 1962-1971, [hep-ph/9707333]. 

[68] G. Belanger, F. Boudjema, A. Pukhov, and A. Semenov, Comput.Phys.Commun. 
180 (2009), 747-767, [0803.2360]. 

[69] R. Kappl, M. Ratz, and M. W. Winkler, Phys.Lett. B695 (2011), 169-173, 
[1010.0553]. 

[70] Super-Kamiokande Collaboration, T. Tanaka et al., Astrophys.J. 742 (2011), 78, 
[1108.3384]. 

[71] R. Kappl and M. W. Winkler, Nucl.Phys. B850 (2011), 505-521, [1104.0679]. 

[72] IceCube Collaboration, R. Abbasi et al, Phys.Rev. D85 (2012), 042002, 
[1112.1840]. 

[73] Fermi-LAT collaboration, M. Ackermann et al., Phys.Rev.Lett. 107 (2011), 241302, 
[1108.3546]. 

[74] PAMELA Collaboration, O. Adriani et al., Phys.Rev.Lett. 105 (2010), 121101, 
[1007.0821]. 

[75] F. Donato, D. Maurin, P. Brun, T. Delahaye, and P. Salati, Phys.Rev.Lett. 102 
(2009), 071301, [0810.5292]. 

[76] K. Abe, H. Fuke, S. Haino, T. Hams, M. Hasegawa, et al., Phys.Rev.Lett. 108 
(2012), 051102, [1107.6000]. 

[77] R. Kappl and M. W. Winkler, Phys.Rev. D85 (2012), 123522, [1110.4376]. 

[78] W. Porod, Comput.Phys.Commun. 153 (2003), 275-315, [hep-ph/0301101]. 

[79] W. Porod and F. Staub, (2011), 1104.1573. 

[80] F. Staub, (2008), 0806.0538. 

[81] F. Staub, Comput.Phys.Commun. 181 (2010), 1077-1086, [0909.2863]. 

20 



[82] F. Staub, Comput.Phys.Commun. 182 (2011), 808-833, [1002.0840]. 
[83] F. Staub, (2012), 1207.0906. 

[84] D. M. Pierce, J. A. Bagger, K. T. Matchev, and R.-j. Zhang, Nucl.Phys. B491 
(1997), 3-67, [hep-ph/9606211]. 

[85] F. Staub, W. Porod, and B. Herrmann, JHEP 1010 (2010), 040, [1007.4049]. 

[86] A. Dedes, G. Degrassi, and P. Slavich, Nucl.Phys. B672 (2003), 144-162, [hep- 
ph/0305127]. 

[87] A. Dedes and P. Slavich, Nucl.Phys. B657 (2003), 333-354, [hep-ph/0212132]. 

[88] A. Brignole, G. Degrassi, P. Slavich, and F. Zwirner, Nucl.Phys. B643 (2002), 
79-92, [hep-ph/0206101]. 

[89] A. Brignole, G. Degrassi, P. Slavich, and F. Zwirner, Nucl.Phys. B631 (2002), 
195-218, [hep-ph/0112177]. 

[90] M. Spira, A. Djouadi, D. Graudenz, and P. Zerwas, Nucl.Phys. B453 (1995), 17-82, 
[hep-ph/9504378]. 

[91] G. Belanger, F. Boudjema, A. Pukhov, and A. Semenov, Comput.Phys.Commun. 

176 (2007), 367-382, [hep-ph/0607059]. 

[92] G. Belanger, F. Boudjema, A. Pukhov, and A. Semenov, Comput.Phys.Commun. 

177 (2007), 894-895. 

[93] G. Belanger, F. Boudjema, A. Pukhov, and A. Semenov, (2010), 1005.4133. 
[94] A. Pukhov, (2004), hep-ph/0412191. 

[95] G. Belanger, N. D. Christensen, A. Pukhov, and A. Semenov, Com- 
put.Phys.Commun. 182 (2011), 763-774, [1008.0181]. 

[96] C. Y. Pages, 'SM Higgs Branching Ratios and Partial-Decay Widths, 
[twiki . cern. ch / twiki /bin / view /LH CPhysics / CERN YellowReport PageBR2] . 

[97] ATLAS Collaboration, ATLAS-CONF-2012-154. 

[98] CMS Collaboration, CMS-SUS-12-022. 



21 



Input 


tan/3 


1.2 


v s [GeV] 


-4.0 


A 


0.74 


A x [GeV] 





K 

lis [GeV] 


1.4 
103.0 


A K [GeV] 
b s [GeV 2 ] 
b/i [GeV 2 ] 
M 2 [GeV] 




3.356 • 10 5 


li [GeV] 
M 1 [GeV] 


280.0 
1500.0 


2.4 • 10 5 
193.0 


M 3 [GeV] 


1500.0 


m S caiar [GeV] 
£s [GeV 3 ] 


1500.0 


A top Y top [GeV] 


1500.0 


0.0 


CP even Higgs sector 


m hl [GeV] 


125.7 


down fraction hi 


41.5% 


m h2 [GeV] 


690.1 


up fraction hi 


57.8% 


m^ 3 [GeV] 


786.8 


singlet fraction hi 


0.7% 




CP odd Hi 


ggs sector 




[GeV] 


247.5 


singlet fraction Ai 


99.9% 


m4 2 [GeV] 


691.9 


up and down fraction Ai 


0.1% 


Neutralino sector 


X°i [GeV] 


130.0 


bino fraction Xi 


<0.1% 


X§ [GeV] 


156.4 


wino fraction Xi 


5.1% 


X§ [GeV] 


316.2 


down-higgsino fraction Xi 


0.3% 


Xl [GeV] 


331.6 


up-higgsino fraction Xi 


10.0% 


$ [GeV] 


1497.4 


singlet fraction Xi 


84.5% 


Chargino sector 


Xt [GeV] 


154.8 


wino fraction xt 


70.6% 


X 2 + [GeV] 


332.6 


higgsino fraction xt 


29.4% 


Electroweak observables 


Rzz 


1.2 
1.0 


Rbb 

Rtt 


1.0 
1.0 


Br(6 — > 57) 


3.4- 10" 4 


By(B s fifi) 


3.7- 10- 9 




-1.2 • lO" 11 


Sp 


4.5 • 10~ 5 


Dark matter 


Vlh 2 


0.1 


Xfo 


24.9 


cr^ 1 [cm 2 ] 
(crv) 77 [cm 3 /s] 


2.2 • 10~ 45 


<7p D [cm 2 ] 
(crv) 7Z [cm 3 /s] 


3.8 • 10~ 40 


0.83 • 10- 27 


0.79 ■ 10~ 27 



Table 1: Benchmark point for the GNMSSM. Rxy denotes the production cross section of 
pp — > h — > XY normalised to the SM expectations (based on the values of the CERN yellow 
pages [96]). 
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